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Abstract: The rare semileptonic B c — > D s ^ £ + £~ decays are studied in the universal extra 
dimension with a single extra dimension scenario. The sensitivity of differential and total 
branching ratios, polarization asymmetries of final state leptons to the compactification 
parameter is presented, both for muon and tau decay channels. Comparing with the 
standard model, the obtained results indicate that there are new contributions to the 
physical observables. Considering the ability of available experiments, it would be useful 
to study these effects. 
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1 Introduction 

The decays induced by flavor changing neutral current (FCNC) b — > s,d transitions play 
an important role to test the standard model (SM) and also they are very sensitive to the 
new physics. The decays of B u ^ s mesons have been intensively studied, and it will make 
the B physics more complete if similar decays of B c meson are included. Considering the 
experimental facilities, it is possible to investigate B c decays with present ability of the 



In the rare B meson decays, new physics contributions appear through the modification 
of the Wilson coefficients existing in the SM or by adding new structures in the SM effective 
Hamiltonian. Among the various extensions of the SM, extra dimensions are specially 
attractive because of including gravity and other interactions, giving hints on the hierarchy 
problem and a connection with string theory. 

The models with universal extra dimensions (UED) allow the SM fields to propagate 
in all available dimensions [1-4]. The extra dimensions are compactified and the compact- 
ification scale allows Kaluza-Klein (KK) partners of the SM fields in the four-dimensional 
theory and also KK excitations without corresponding SM partners. Throughout the UED, 
a model including only a single universal extra dimesion is the Appelquist-Cheng-Dobrescu 
(ACD) model [5]. The only additional free parameter with respect to the SM is the inverse 
of the compactification radius, 1/R. In particle spectrum of the ACD model, there are 
infinite towers of KK modes and the ordinary SM particles are presented in the zero mode. 

Discussion on the free parameter, both experimental and theoretical aspects of the 
subject, shows that the lower bound should be between 1/R> 250 GeV or 1/R > 350 GeV 
[6-11]. Therefore, we will pay more attention to 1/R = 250 — 350 GeV region but also 
consider 1 / R = 200 — 1000 GeV in our calculations. 

The effective Hamiltonian of several FCNC processes [12]-[13], semileptonic and radia- 
tive decays of B mesons [14-22] and FCNC baryonic decays [10, 23-25] have been inves- 
tigated in the ACD model. Polarization properties of final state particles in semileptonic 
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decays, which is an powerful tool in searching new physics beyond the SM, have also been 
studied widely besides the other observables in these works, e.g [11, 14, 18, 22]. 

The main aim of this paper is to find the possible effects of the ACD model on some 
physical observables related to the B c — > (D s , D)£ + £~ decays. We study differential decay 
rate, branching ratio, and polarization of final state leptons, including resonance contribu- 
tions in as many as possible cases. We analyze these observables in terms of the compacti- 
fication factor and the form factors. The form factors for B c — > (D s , Z?)£ + ^~processes have 
been calculated using different quark models [26-29] and three-point QCD sum rules [30] . 
In this work, we will use the form factors calculated in the constituent quark model [26]. 

This paper is organized as follows. In section 2, we give the effective Hamiltonian for 
the quark level processes b — > s,d£ + £~ with a brief discussion on the Wilson coefficients in 
the ACD model. We drive matrix element using the form factors and calculate the decay 
rate in section 3. In section 4, lepton polarizations are evaluated and in the last section, 
we introduce our numerical analysis and discuss the obtained results. 



2 Theoretical Framework 

The effective Hamiltonian describing the quark level b — > (s,d)£ + £~ processes in the SM 
is given by [31] 



eff 



y/2ir 



v tq >v t t 



C e g ff (q'^L b) hn + C w (q'^L b) £>y» l5 £ 



2C e 7 Sf m h {q'ia^Rb)hn 



(2.1) 



where q = Pb c — PD q , is the momentum transfer and q' = s,d. 

New physics effects in the ACD model come out by the modification of the SM Wilson 
coefficients appear in the above Hamiltonian. This process can be done by writing the 
Wilson coefficients in terms of 1/R dependent periodic functions the details of which can 
be found in [12, 13]. That is, Fo(xt) in the SM is generalized by F(xt, l/R) accordingly 



F(x t , l/R) = F (xt) + J2 x ^ 



(2.2) 



n=l 



where xt = m^/m^, x n = m^/m^, and m n = n/R. These modified Wilson coefficients 
can widely be found in the literature. Here, we will not discuss the details and follow the 
explanation given in [11]. 

In the ACD model, a normalization scheme independent effective coefficient can 
be written as 

C e 7 ff (fi b ,l/R) = r, 1 V 23 C 7 (fi w ,l/R) 



+ t:(v U/23 ~ V 16/23 )C 8 (W, l/R) + C 2 (nw, l/R) E • (2-3) 



-2- 



The coefficient has perturbative part and we will also consider the resonance 
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contributions coming from the conversion of the real cc into lepton pair. So, C% is given 



by 

C 9 e/ V, l/R) = C 9 ( M , 1/12) (l + ^w(s')) + Y(fj) + C 9 res (m) • (2.4) 

For Cg, in the ACD model and in the naive dimensional regularization (NDR) scheme we 
have 



Cgfji, l/R) = P NDR + Y{Xt, 2 \ /R) - 4Z(x t , l/R) + P E E(x t , l/R) (2.5) 



where Pq DR = 2.60 ±0.25 and the last term, PgEixt, l/R), is numerically negligible. The 
perturbative part, coming from one-loop matrix elements of the four quark operators, is 

Y(ji) = %,a)[3Ci(At) + C 2 (ji) + 3C 3 (m) + C7 4 (m) + 3C 6 (aO + C 6 (aO] 
- ^(1, a) (4C 3 (aO + 4C 4 (/i) + 3C 5 (aO + C 6 (/i)) 
-\h(0, s ) [<7 3 ( M ) + 3C 4 ( M )] 

+ H (3C 3 (/i) + C 4 (m) + 3C 5 (^) + C 6 (m)) , (2.6) 
y 

with y = m c /mb. The explicit forms of the functions appear above equations can be found 
in [32, 33]. The resonance contribution can be done by using a Breit-Wigner formula [34] 

_ 3 ^ nm z> ^T>^ 

x [3Ci(/i) + C7 2 (m) + 3C 3 (^) + C 4 (/i) + 3C 5 (aO + CeM] ■ (2.7) 

The normalization is fixed by the data in [35] and k is taken 2.3. 

The Wilson coefficient C\q is independent of scale fi and given by 

_ iw/g) 

WO • (2.SJ 

3 Matrix Elements and Decay Rate 

The matrix elements for B c — > D g /£ + £~ can be written in terms of the invariant form 
factors over B c and D q >. The parts of transition currents containing 75 do not contribute, 
so the non- vanishing matrix elements are 



Mq 2 



Dq'{PD gl ) \q'icr^q v b\ B c (p Bc )) = ~ ( p B c + Pog'W - q^( m B, 



111 



D 



D q >( P D,) \q%b\ B c (p Bc )) = f + (q z ){P Bc + Pn ,)„ + f-tfK ■ (3.1) 
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Figure 1. (color online) The dependence of differential branching ratio on s without resonance 
contributions for B c -> D s l + i~. (In the legend l/R = 200, 350, 500 GeV.) 




Figure 2. (color online) The dependence of differential branching ratio on s without resonance 
contributions for B c — > Di + l~ . 



The transition amplitude of the B c — > D q >£ + £ decays can be written using the effective 
Hamiltonian and eq. (3.1) as 



M {B c -> D q 



A(P Bc +P D ,)„+Bq, 



+ hV C{P Bc +P Dal )^ + D qtl 



(3.2) 



with 



A = C e 9 ff f + + 



2m b f>. 



B = C e 9 ff f. 



C — C10/+ 
D = C 10 f- 



m Bc + m D , 



T n eff 



2m b (m 2 Bc -m%j)fr 
q 2 (m Bc + m D ,) 



-a 



eff 
7 ' 



(3.3) 
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Figure 3. (color online) The dependence of differential branching ratio on s with resonance 
contributions for B r — > D*i + i~ . 
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Figure 4. (color online) The dependence of differential branching ratio on s with resonance 
contributions for B r D£ + £~ . 



Finally, following dilepton mass spectrum is obtained by eliminating angular depen- 
dence in the double differential decay rate, 

dT G 2 a 2 m Bc , .2 „ /T a fi a\ 

where s = q 2 /m 2 Bc , A = 1 + r 2 + s 2 — 2r — 2s — 2rs, r = m 2 D , /m 2 B ^, v = ■Jl — Am 2 / sm 2 Bc 
and 
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A Dql = -m Bc (3 - v 2 )X(\A\ 2 + \C\ 2 ) + 4m% c s(2 + r-s)(l- v 2 ) |Cf 
+ 16m| mjs \D\ 2 + 32m| mj(l - r)Re(CD*) . 



(3-5) 



4 Lepton Polarization Asymmetries 

Studying polarization asymmetries of final state leptons is an useful way of searching new 
physics. Therefore, we will discuss the possible effects of the ACD model in the lepton 
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Figure 5. (color online) The dependence of branching ratio on 1/R with and without resonance 
contributions for B c — > D s £ + £~ . (The subscript R represents resonance contribution.) 
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Figure 6. (color online) The dependence of branching ratio on 1/R with and without resonance 
contributions for B r — > D£ + £~. 



polarization. Using the convention in [37, 38], we define the orthogonal unit vectors S~ 
(i = L, T, N) for the longitudinal, transverse and normal polarizations in the rest frame of 
l~ as 



Sl = (0,e L ) 



0. 



Pi 

m 



S T = (0, e T ) = (0, e N x e L ) , 



7 N 



(0,e N ) 



PD ql x Pe 
'PD q , x Pi 



(4.1) 



where pe and po q , are the three momenta of £~ and D q > meson in the center of mass (CM) 
frame of final state leptons, respectively. The longitudinal unit vector 5^ is boosted by 
Lorentz transformation, 



q -» _ i'Wi\ EePe 

L, CM ~ 1 



m e mi\pi\ 

while vectors of perpendicular directions remain unchanged under the Lorentz boost. 



(4.2) 
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Figure 7. (color online) The dependence of longitudinal polarization on s without resonance 
contributions for B c — > D s £ + £~. 
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Figure 8. (color online) The dependence of longitudinal polarization on s without resonance 
contributions for B c — > D£ + £~. 



The differential decay rate of B c — > D q i 
be written in the following form 



for any spin direction n of the i can 



dT(ri-) _ 1 (dT\ 



ds 



2 \ds J. 



1+ I P L e L +P AT e M + Pr r e. 



N N ~ 1 T^T 



■ n 



(4.3) 



where (dT/ds)o corresponds to the unpolarized decay rate, the explicit form of which is 
given in eq. (3.5). 

The polarizations P£, P^ and PZ in eq. (4.3) are defined by the equation 



pr(s) 



dT(n~ = er)/ds - dT(n~ = -e,r)/ds 



dT(n- = er)/ds + dT{n 



-e, )/ds 



Here, PT and P^ represent the longitudinal and transversal asymmetries, respectively, of 
the charged lepton £~ in the decay plane, and PZ is the normal component to both of 
them. 
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Figure 10. (color online) The dependence of longitudinal polarization on s with resonance 
contributions for B c — > Di + l~ . 



The straightforward calculations yield the explicit form of the longitudinal polarization 
for B c ^ D sd £+£- as 



1 6 



(4.4) 



and the transversal polarization is given by 



Pr, 



A 



(r-1) 



Re[AC*] + Re[AD*\ 



(4.5) 



The normal component of polarization is zero so we have not stated its explicit form 
here. 



5 Numerical Analysis and Discussion 

In this section, we will introduce numerical analysis of physical observables. The input 
parameters used in this work are me c = 6.28 GeV, m£, a = 1.968 GeV, = 1.870 GeV, 
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Figure 11. (color online) The dependence of longitudinal polarization on 1/R with resonance 
contributions for B r — > D*£ + £~ . 
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Figure 12. (color online) The dependence of longitudinal polarization on 1/R with resonance 
contributions for B r —t Dl + £~ . 



m b = 4.8 GeV, m M = 0.105GeV, m T = 1.77 GeV, \V tb V t * s \ = 0.041, \V tb V* d \ = 0.008, 
G F = 1.17 x 10" 5 GeV' 2 and t Bc = 0.46 x 10~ 12 s [35]. 

To make numerical predictions, we also need the explicit forms of the form factors 
/+,/_ and fa. In our analysis, we used the results of [26], calculated in the constituent 
quark model and q 2 parametrization is given by 



vy > i- a ( q ym 2 Bc ) + b( q ym 2 Bc y ' 

where the values of parameters -F(O), a and b for the B c — > (D S ,D) decays are listed in 
Table 1. 

In this work, we also took the long-distance contributions into account. While doing 
this, to minimize the hadronic uncertainties we introduce some cuts around J /ip and ip(2s) 
resonances as discussed in [11]. 

In the analysis, first the differential branching ratios are calculated with and without 
resonance contributions and s dependence for 1/R = 200, 350, 500 GeV are presented in 
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Figure 13. (color online) The dependence of transversal polarization on s without resonance 
contributions for B c — > D s £ + £~. 




Figure 14. (color online) The dependence of transversal polarization on s without resonance 
contributions for B c — > D£ + £~ . 

figures 1-4 for B c — > (D s , D)£ + £~ . One can notice, the change in differential decay rate and 
difference between the SM results and new effects in the figures. The maximum deviation 
is around s = 0.32(0.39) in figure 1 and s = 0.36(0.40) in figure 2 for /i(r). At these s 
values, the deviation is ~ 50% more than that of the SM results for 1/R = 200 GeV in 
all decay channels and ~ 20% for 1/R = 350 GeV. For 1/R > 500 GeV, the deviation 
becomes ~ 10% and less. Considering the resonance effects, the differential decay rates 
also differ from their SM values as 1/R — > 200 GeV, which can be seen in figures 3 and 4. 
So, studying differential decay rate, particularly in 1/R = 250 — 350 GeV region, can be 
an appropriate tool for searching the effect of extra dimension. 

To introduce the contributions of the ACD model on the branching ratio, we present 
1/R dependent ratios with and without resonance cases in figures 5 and 6. The common 
feature is that as 1/R increases, the branching ratios approach to their SM values. In all 
decay channels, for 1/R ~ 500 GeV the deviations are more than 10% from their SM values. 
Additionally, around 1/R ~ 200 GeV the ACD contribution is more than 50%. Keeping the 
discussion on lower bound on 1/R in mind once again, in the region 1/R ~ 350 — 250 GeV 
the branching ratios increase with an average of (25 — 35)%. So, the calculated branching 
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Figure 16. (color online) The dependence of transversal polarization on s with resonance contri- 
butions for B c -> D£ + £~. 



ratios without resonance contributions for the SM and in between these bounds we find, 

(1.24,1.47 - 1.67) x 1(T 7 
(0.29,0.36 -0.41) x 10~ 7 
(0.29,0.34 - 0.39) x 10~ 8 
(0.077,0.092 - 0.106) x 10~ 8 . 

Here, the first value in any branching ratios above is corresponding to the SM, while 
the rests are for 1/R = 350 — 250 GeV . A similar behavior is valid for resonance case which 
can be followed by the figures. 

Adding the uncertainty on the form factors may influence the contribution range of 
the ACD model. However, the variation of the branching ratio, calculated with the central 
values of form factors, in the ACD model for different 1/Rs with the SM values, can be 
considered as a signal of new physics. 

The dependence of longitudinal polarization on s without resonance contributions are 
given by figures 7-8. As 1/R approaches to 200 GeV, the longitudinal polarization differ 



Br{B c -> D sI i + ll-) = 
Br(B c -> D s t+t-) = 
Br(B c -> D fi+fT) = 
Br(B c -> D t + t~) = 
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Figure 17. (color online) The dependence of transversal polarization on 1/R with resonance 
contributions for B r — > D*£ + £~ . 
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Figure 18. (color online) The dependence of transversal polarization on 1/R with resonance 
contributions for B c — > D£ + £~ . 



from the SM values, slightly. This is more clear in 0.34 (0.35) < s < 0.45 (0.47) for 
B c — > D s (D)t + t~ decay, respectively. For // channels, the deviation is valid in all s 
range, excluded the minimal and maximal points of s. In B c — > D s (D)fi + decay, the 
difference is, for example at s ~ 0.38, 5% (3%), respectively, for 1/R = 200 GeV. For 
B c — > D s (D)t + t~ decay the results are more significant that is at s ~ 0.40(0.42), the 
variation for 1/R = 200 GeV is 12% (8%), respectively. Including resonance contributions, 
in addition to the above effects, deviation between the resonant and nonresonant values 
increases as l/R—> 200 GeV (Figures 9 and 10). 

In order to clarify the dependence on 1/R, we eliminate the dependence of the lepton 
polarizations on s, by considering the averaged forms over the allowed kinematical region 
as 



{Pi) 



(l-m D(j ,/m Sc ) 2 fig 

Pi—ds 

{2 m t /m Bc )l as 

(l-m D(j/ /m Bc ) 2 d g 

—ds 

(2m,/m Bc ) 2 ds 



(5.1) 
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The 1/R dependant average longitudinal polarizations are given in figures 11 and 12. As 
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Table 1. B c — > D S ( j decays form factors calculated in the constitute quark model. 

it can be seen from the figures, the maximum deviation is 2% for fi channels and 9% (6%) 
for B c -> D s (D)t + t~, respectively, at 1/22 = 200 GeV. 

The variation of transversal polarization with respect to s are given by figures 13-16. 
In fj, channels the difference is negligible whereas in r channels up to s ~ 0.46 (0.48) for 
B c —7- D s {D)i + i~ decays, respectively, the effects of the UED can be seen. At, for example, 
s ~ 0.40(0.44), the SM value vary 5 - 12% (4-9%) for 1/R = 350-200GeV, respectively. 
Finally, the average transversal polarization can be followed by figures 17 and 18. For the 
\x channels the maximum deviation is 2% and for B c — > D s (D)t + t~ that is 10%(6%) at 
1/R = 200 GeV. 

As an overall conclusion, as 1/R — > 200 GeV the physical values differ from the SM 
results. Up to a few hundreds GeV above 250 GeV or 350 GeV, which are more common 
lower bounds, one can see the effects of UED. Considering the differential and integrated 
branching ratios, there appears essential difference, particularly in between 1/R = 250 — 
350 GeV comparing with the SM results. The polarization effects in \i channels are either 
irrelevant nor negligible while in r channels small effects are obtained. Under the discussion 
in this work, studying these decays experimentally can be useful. 
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